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Small lanthanide clusters have interesting magnetic properties, but their structures are unknown. We have 
identified the structures of small terbium cluster cations Tb+ (n = 5 — 9) in the gas phase, by analysis of their 
vibrational spectra. The spectra have been measured via IR multiple photon dissociation of their complexes 
with Ar atoms in the 50 — 250 cm -1 range with an infrared free electron laser. Density functional theory 
calculations using a 4/-in-core effective core potential (ECP) accurately reproduce the experimental far-IR 
spectra. The ECP corresponds to a 4/ 8 5<i 1 6s 2 trivalent configuration of terbium. The assigned structures are 
similar to those observed in several other transition metal systems. From this, we conclude that the bonding 
in Tb clusters is through the interactions between the bd and 6s electrons, and that the 4/ electrons have 
only an indirect effect on the cluster structures. 



Small clusters of lanthanide elements are promising ma- 
terials for molecular magnets, due to their large mag- 
netic moments, and strong spin-orbit interaction.^ The 
magnetism of lanthanide metal clusters originates from 
the indirect (RKKY) exchange interaction, and the mag- 
netic anisotropy is due to the interaction between the 
asymmetric 4/ charge cloud and the local crystal field. 
Both of these properties are highly dependent on the geo- 
metric structure, and the lack of knowledge of the cluster 
structures is a major barrier to understanding their mag- 
netism. To obtain fundamental insights into their prop- 
erties, the clusters can be studied in gas-phase molecular 
beam experiments, which provide an environment free of 
external influences. 

The magnetism of lanthanide clusters has been stud- 
ied experimentally using the Stern-Gerlach deflection 
technique. 2 These experiments demonstrated that lan- 
thanide clusters 3 have large size-dependent magnetic 
moments and in some cases, Curie temperatures much 
higher than the corresponding bulk materials. Theoreti- 
cally, the Gdi3 cluster has been most extensively studied: 
In Ref. HJ the cluster structures were taken to be frag- 
ments of the bulk hep lattice, and the magnetism was 
described phenomenologically with a Heisenberg Hamil- 
tonian. The exchange constants which describe the 
strength of the interaction between local 4/ moments due 
to the RKKY interaction were taken as undetermined pa- 
rameters, rather than properly estimated from electronic 
structure calculations. In fact, little to nothing is known 
about the physics of the RKKY interaction in small sub- 
nanometer particles. This situation stands in sharp con- 
trast to the present understanding of the bulk lanthanides 
where many physical properties are well understood with 
simple models, and first-principles theory can calculate 



magnetic ordering temperatures with near quantitative 
accuracy.^ The structures of the bulk lanthanides are 
all known from X-ray diffraction on single crystal and 
thin film samples, and the lack of such structural infor- 
mation is a major barrier preventing the application of 
high-level theory to small clusters. A recent theoretical 
study of the Gdi3 cluster 7 which attempted to include the 
4/ electrons concluded that the lowest energy structures 
were icosahedral, rather than /icp-like, contradicting the 
previous results of. Refs. |4] Despite this, lanthanide clus- 
ters have received little or no attention from experimental 
structural studies. 

State of the art computational electronic-structure the- 
ory still has many difficulties with reliable prediction of 
cluster structures and their related properties from first 
principles. This holds in particular for metals with par- 
tially filled d (transition metals) or / (lanthanides, ac- 
tinides) shells or for heavy elements, like gold, where 
relativistic effects are significant. Experimental data is 
therefore essential for testing the validity of predictions 
and in motivating methodological developments. To this 
end, hybrid experimental/theoretical studies with photo- 
electron spectroscopy^ trapped- ion electron diffraction,^ 
and far-IR spectroscopy^"^ have provided invaluable in- 
formation on the structures of small metal clusters. To 
identify the cluster structures, we present here exper- 
imental far-IR vibrational spectra for cationic terbium 
clusters in the gas phase, mass-selectively measured by 
dissociation of weakly bound Ar messenger complexes. 
This method has been suc cessfully ap plied before, e.g., 
to transition metal clusters!- 10 * 11 * 1 SH UT Terbium was cho- 
sen because there is a single naturally occurring isotope, 
which is beneficial for mass spectrometric studies, and it 
can be reasonably expected to behave as a prototype for 
the heavy lanthanides 5 (Gd, Tb, Dy, Ho, Er), which are 
all trivalent and have the same hep bulk lattice type. 
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The Free Electron Laser for Infrared experiments (FE- 
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LIX) in Nieuwegein, the Netherlands, was used as the 
far-IR source for these experiments. The experimen- 
tal setup a nd me thods have been described in previous 
publication d 10 ! 13 ! and we only give a brief summary of 
the relevant details here. Tb clusters are produced in a 
laser vaporization cluster source which is in part cooled 
to 77 K. The He carrier gas is seeded with 5 % Ar, which 
leads to the formation of a small population of weakly 
bound Ar complexes of the Tb+ clusters. The cluster 
beam is irradiated by the tunable far-IR radiation from 
FELIX and mass analyzed with a reflectron time-of-flight 
mass spectrometer. The FELIX beam used in the exper- 
iment delivers an energy on the order of ~ 10 mJ within 
a pulse of about 7-9 /is duration. Mass spectra of alter- 
nating shots from the cluster source with FELIX toggled 
on and off are stored into separate channels of a digitizer 
and in total about 500 mass spectra are averaged per 
channel and frequency interval. The IR induced relative 
depletion of the mass signals for the Tb+-Ar complexes 
is evaluated as a function of IR frequency to obtain the 
size specific far-IR spectra. 

The experimental IR spectra for Tb+ (n = 5 — 9) are 
shown in Figs. [T]- [3] There are well-resolved absorption 
bands for most cluster sizes. The spectra of the smaller 
sizes, n = 5, 6, and 7, are dominated by a single intense 
feature, which may suggest a highly symmetric struc- 
ture, while the larger sizes, n = 8 and 9 have more richly 
featured spectra. For the studied cluster sizes we find 
absorption bands in the 50-160 cm -1 range, which is ap- 
proximately the range where the optical phonon branches 
of the bulk metal were observed by neutron scattering.^ 
With few exceptions, all observed metal cluster vibra- 
tions are below the reported value for the stretching fre- 
quency of the neutral Tb2 of 137 cm _1 P^ 

To identify the clusters' structures based on the exper- 
imental spectra we perform quantum chemical calcula- 
tions using density functional theory (DFT). We use 
the Perdew-Burke-Ernzerhof (PBE) functional for ex- 
change and correlation as implemented in Turbomole 
6.20 The treatment of the 4/ electrons with DFT re- 
mains a challenge because they occupy an open shell and 
are highly localized to their parent atoms, where corre- 
lation effects are strong. Their strong localization also 
implies that they do not directly participate in chemi- 
cal bonding. Standard DFT functionals are known to 
falsely delocalize the 4/ electrons, and explicit treatment 
of the 4/ states requires corrections such as LDAH-UP^ 
Alternatively, it is possible to exploit the localization of 
the 4/ states and treat them as a part of the core with 
an effective core potential. We follow this approach, us- 
ing the ECP54MWB effective core potential developed 
by Dolg et al. 22 The valence orbit als are described with 
the (6s6p5d)/[4s4p4d] + 2slpld (ECP54MWB-II) ba- 
sis set 23 . The ECP includes 8 4/ electrons in the core. 
This corresponds to the commonly observed trivalent 
(4/ 8 5<i 1 65 2 ) configuration of the terbium atom. 

For metal clusters the interrogation of the configurational 
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FIG. 1. Experimental far-IR spectra of Tbjj" and Tbjt" com- 
pared with the theoretically predicted spectra for a selec- 
tion of low-lying structural isomers and 5d6s spin states. 
The experimental spectra are shown in relative cross sections 
a oc v/P(y) log(yj^y) derived from the depletion of the Ar 
complexes, normalized by the photon fluence. All experimen- 
tal spectra are obtained from the complexes with a single Ar 
atom, except for Tb^~ where for intensity reasons the complex 
with two Ar atoms has been used. To help the eye, the exper- 
imental raw data points (gray dots) are overlayed with their 
5 point running average. For each isomer its energy relative 
to the lowest energy isomer identified (in eV) and the num- 
ber of unpaired electrons 2S are given. The calculated stick 
spectrum is folded with a Gaussian line width function of 2 
cm -1 full width at half maximum. 



space becomes rapidly more complex with growing size 
and many locally stable arrangements of atoms, in addi- 
tion to the true global minimum, may be identified for 
any given cluster size. We investigated a wide range of 
possible isomers starting our search with geometries that 
have been reported for transition metal clusters (after 
scaling the interatomic distances to values which match 
the bulk values for Tb) including those tha t have been 
identified in our previous studies JMnHHHl followed by 
local optimization. In addition, for several sizes, we used 
Monte-Carlo basin-hopping sampling 24 to locate other 
local minima. No symmetry constraints have been ap- 
plied. For all local minima found the IR spectra have 
been calculated analytically based on the harmonic ap- 
proximation; the frequencies are not further scaled. For 
n = 5 — 7, alternative isomers are typically much higher 
in energy. This suggests that the potential-energy surface 
is quite smooth. 

An additional parameter of the calculations is the elec- 
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FIG. 2. Far-IR spectrum and comparison to predictions for 
different isomers of Tb^ . Other plausible structures for Tb^ , 
such as a capped octahedron, or bicapped pentagonal bipyra- 
mid are either significantly higher in energy, or converge to 
a pentagonal bipyramid under local optimization of the total 
energy. 



tronic occupation of the valence states. In all of these 
calculations electronic states with different numbers of 
unpaired 5d/6s valence electrons have been considered 
in order to identify the preferred configuration. The 
spin states we refer to in the following only relate to 
the bd/6s occupation as the 4/ electrons are part of the 
ECP and not explicitly treated. If the 4/ electrons were 
included, the total spin of any Tb cluster would be much 
higher. (The 4/ shell of trivalent Tb has a configura- 
tion J = 6,gL = 3/2 yielding a theoretical moment of 
Ql^bJ = 9/iB per atom!^) The ECP calculations also ex- 
plicitly neglect any ordering of the 4/ moments, as well as 
the 4/ — (5d6s) exchange interaction, which is known to 
split the majority and minority (5d6s) spin states for bulk 
Tb metal in the ferromagnetic state by 0.1 eV! 2 ^ An ex- 
change splitting of this magnitude could plausibly change 
the ground state valence occupation, and thus the Far-IR 
spectra. The surprising fact is that the 4/-in-core model, 
which completely neglects this effect, still gives remark- 
ably good agreement with the experimental spectrum. 
The frequencies of the IR active bands are correctly pre- 
dicted to within 1%, and the relative intensities match 
qualitatively. Figs. [I|j3] show a comparison between the 
experimental far-IR spectra and those of the identified 
putative ground states. More isomers and their spectra 
calculated with a smaller, presumably less accurate, ba- 
sis set are reported in the supporting information 2 ^. In 
short, we can assign the following structures for Tb^lg 
(approximate symmetries are given in parentheses): 5 - 
trigonal bipyramid {C^v 2*5 = 4), 6 - distorted octahe- 
dron (6c - D2h 2S = 1), 7 - pentagonal bipyramid or 
decahedron (D 5 h 2£ = 2), 8 - bicapped octahedron (D 2 d 
2S = 3), 9 - tricapped octahedron (C2 V 25 = 2). 

For the small clusters, Tb§_ 7 , the best matches to the 
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FIG. 3. Far-IR spectrum and comparison to predictions for 
different isomers of Tb+ and Tb+. 



IR spectra are obtained with the putative ground states 
which are calculated to have 4, 1, and 2 unpaired 5d/6s 
electrons, respectively. For a given structure, different 
spin states are often close in energy, but their calculated 
IR spectra can differ remarkably, as seen for TbJ. In this 
case, the doublet state 8b is higher in energy than the 
quartet 8a, but the spectrum of 8a fits the experiment 
better. The presence of (several) unpaired 5d/6s elec- 
trons is consistent with Stern- Gerlach experiments on La 
and Y clusters, 26 which found many sizes with net mag- 
netic moments ( 0.2 /i^/atom). La and Y have very simi- 
lar valence electronic structure to the lanthanides, but no 
4/ electrons. Thus the observation of high spin in these 
clusters supports our finding, even without properly ac- 
counting for the possibility of a 5d6s spin polarization 
induced by the large localized 4/ moments. 

We find no evidence of multiple isomers for any of these 
small Tb clusters, although for Tbg" there are several iso- 
mers based on two structural motifs (9a/9b and 9c/9d) 
which are close in energy (within 0.2 eV) and which have 
similar spectra. 9a/ 9b is formed by capping three adja- 
cent faces of an octahedron and 9c/9d can be described 
as two interpenetrating decahedra. Both of these could 
be present, but structure 9a/9b fits the experimental 
spectra better, and the energy difference of 0.1 eV is too 
small to be considered meaningful. 

The lowest energy isomer identified by theory does not al- 
ways give the best match to the theory, a case illustrated 
by Tb^j". Tb^j" has 17 (5d6s) electrons, one fewer than 
the 18 needed to satisfy a spherical shell closing. This 
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would indicate a preference for a low-spin distorted ge- 
ometry 6c, over the octahedral 6a. However, the highly 
symmetric octahedral structure has high density of states 
at the HOMO, which increases the energy gained by ex- 
change splitting in a high spin state. It is possible for this 
effect to stabilize a highly symmetric structure against a 
Jahn- Teller distortion. Indeed, the calculations predict 
that the 2S = 5 octahedral structure 6a is 0.3 eV lower 
in energy than the distorted low-spin (2S = 1) structure 
6c, but 6c appears to match the experimental spectrum 
better. It is important to note, however, that the geo- 
metric difference between 6a and 6c is small, (the bond 
lengths differ by 2-5%) and their vibrational spectra are 
very similar. The octahedral structure 6a actually has a 
triply degenerate vibrational mode located at 100 cm -1 , 
which is close to where an intense feature is observed in 
the experimental spectrum. However, the IR intensities 
calculated for these modes are very low. Further investi- 
gation is required to understand whether the presence of 
the Ar, or the inclusion of the 4/ electrons would change 
the IR intensities of these modes. Thus our experiment 
cannot determine conclusively whether Tbg" is a perfect 
octahedron. 

To conclude, we have demonstrated that far-IR spec- 
troscopy allows for structural assignment for small lan- 
thanide clusters. Our theoretical approach, where the 4/ 
electrons are treated as a part of the core, predicts the 
frequencies and IR intensities of the cluster vibrational 
modes with remarkable accuracy. This confirms that, 
like the bulk materials, the 4/ electrons remain tightly 
localized in metal clusters and do not directly partici- 
pate in bonding. The structures identified are similar 
to the geometries found for transition metal clusters. In 
this sense, heavy lanthanides behave like early transition 
metals with small d-occupation. It is also highly sur- 
prising that the structures can be determined while the 
magnetic order and 4/ — (5d6s) exchange interaction are 
completely neglected. Knowledge of the cluster struc- 
tures will enable the application of higher level theory 
for an explicit treatment of the 4/ electrons, which in 
turn will provide a detailed understanding of lanthanide 
cluster magnetism. 
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